Cluster cannibalism and scaling relations of galactic stellar nuclei 
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(M ■ ABSTRACT 

Q . Recently, very massive compact stellar systems have been discovered in the 

"^ . intracluster regions of galaxy clusters and in the nuclear regions of late-type disk 

■~>^ . galaxies. It is unclear how these compact stellar systems - known as "ultra- 

p . compact dwarf" (UCD) galaxies or "nuclear clusters" (NCs) - form and evolve. 

Q . By adopting a formation scenario where these stellar systems are the product of 

^ . multiple merging of star clusters in the central regions of galaxies, we investigate, 

c^ . numerically, their physical properties. We find that physical correlations between 

velocity dispersion, luminosity, effective radius, and average surface brightness in 
\i/ , the stellar merger remnants are quite different from those observed in globular 

clusters. We also find that the remnants have triaxial shapes with or without 
figure rotation, and these shapes and their kinematics depend strongly on the 
initial number and distribution of the progenitor clusters. These specific pre- 
dictions can be compared with the corresponding results of ongoing and future 
observations of UCDs and NCs, thereby providing a better understanding of the 
origin of these enigmatic objects. 

Subject headings: galaxies: nuclei — galaxies: dwarfs — globular clusters:general 
- galaxies: star clusters 
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Introduction 



A new type of sub-luminous and extremely compact "dwarf galaxy" has been recently 
discovered in an 'all-object' spectroscopic survey centered on the Fornax Cluster (Drinkwater 
et al. 2000). These "dwarf galaxies" , which are members of the Fornax Cluster, have intrinsic 
sizes of less than 100 pc and absolute B band magnitude ranging from —13 to —11 mag and 
are thus called "ultra-compact dwarf" (UCD) galaxies. Although these UCDs are suggested 
to originate from stellar nuclei of bright nucleated dwarf galaxies (Drinkwater 2000; Bekki 
et al. 2001, 2003), it is unclear how such massive nuclei are formed in the central region of 
dwarf galaxies. 

Recent Hubble Space Telescope (HST) photometric observations have discovered very 
luminous nuclear clusters (NCs), with /-band absolute magnitudes (Mj) ranging from —8 
to —14 mag, in the central regions of late- type spirals (Phillips et al. 1996; Carollo et al. 
1998; Matthews et al. 1999; Boker et al. 2002, 2004a,b). The observation that some of the 
bright NCs are quite massive - that is their luminosity does not derive from a small number 
of hot young stars - has raised the question as to how such massive NCs can be formed in 
the central regions of late-type spiral galaxies (Boker et al. 2004a,b). 

One formation scenario for these very massive star clusters (hereafter referred to as 
"VMSCs) is that ordinary star clusters, which can quickly spiral into the nuclear regions 
of galaxies due to dynamical friction, merge with one another to form a single VMSC (i.e., 
galactic stellar nuclei; Tremaine et al. 1975). The physical properties of VMSCs formed 
in this way, however, have not been theoretically/numerically investigated extensively (e.g., 
Fellhauer & Kroupa 2002). In particular, theoretical predictions of the correlations between 
their structural and kinematical properties (e.g., central velocity dispersion) are generally 
considered to be important, because such dynamical correlations (or "scaling relations") 
for a self-gravitating system are generally considered to help discriminate between different 
formation mechanisms (e.g., Djorgovski 1993). In the light of recent discoveries of UCDs and 
NCs, it is thus timely and important to discuss whether such a merger scenario (referred to 
as "cluster cannibalism" in galactic nuclei) is consistent with their observed properties. 

The purpose of this Letter is to provide the first theoretical predictions on the structural 
and kinematical properties of VMSCs formed by dissipationless multiple cluster merging 
based on self-consistent numerical simulations of nucleus formation. We focus particularly on 
correlations between properties such as luminosity (L), effective radius (-Re), central velocity 
dispersion (o-q), and surface brightness at R^. The predicted scaling relations combined with 
current and future observations of UCDs (e.g., Drinkwater et al. 2003) and NCs (e.g., Boker 
et al. 2004a,b) can provide new insight into the origin of galactic stellar nuclei. Dissipative 
formation, which is an alternative formation scenario for VMSCs (Boker et al. 2004a,b), will 
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be discussed in forthcoming papers. 



2. The Model 

We investigate the dynamical evolution of a self-gravitating system composed of smaller 
star clusters (SCs), via numerical simulations carried out on a GRAPE board (Sugimoto et 
al. 1990). Each of the individual SCs that merge with one another to form a VMSC, are 
assumed to have a Plummer density profile (e.g., Binney & Tremaine 1987) with luminosities 
(Lsc) and central velocity dispersions (o"sc) consistent with the relation observed for GCs 
(Djorgovski et al. 1997): 

Lsc oc dsc^''^. (1) 

The scale length (ogc) of a SC is determined by the formula 

a,e = GMjQaJ, (2) 

where G and M^c are the gravitational constant and the mass of the SC, respectively. Since 
the-mass-to-light-ratio {Msc/Lsc) is assumed to be constant for all SCs, Osc and o"sc are 
determined by the equations (1) and (2) for a given Lsc (or Mgc). The normalization factor 
in equation (1) is determined by using the observed typical mass (6 x 10^ ^o), the half 
mass radius (lOpc), and the central velocity dispersion (7 km s~^) for GCs (e.g., Binney & 
Tremaine 1987). 

A SC system is assumed to have either a uniform disk distribution (referred to as "2D" 
models) or a uniform spherical distribution ("3D"), and to be fully self-gravitating (i.e., not 
influenced dynamically by its host galaxy). For the 2D models, a SC system is assumed to 
have only rotation initially (no velocity dispersion), with its rotational velocity (V"rot) at a 
distance (r) from its centre given by: 

V,(r)=CvKir(r), (3) 

where Vcir{r) is the circular velocity at r for this system and Cy is a parameter that de- 
termines how far the system deviates from virial equilibrium (i.e., < Cy < 1). For 3D 
models, a SC system is supported only by its random motions; its one-dimensional isotropic 
dispersion (ctj) at radius r is given by: 



<^^ir)=Cy^-^, (4) 

where U{r) is the gravitational potential of the system. Here a SC system with Cy = 1 
means that the system is in virial equilibrium (or supported fully by rotation). We present 
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the results of models in which all of the SCs are within ~ 100 pc of each other, because they 
are the most consistent with observations of the structural properties of VMSCs. 

We investigate two different representative cases: (1) the SC system is composed only of 
equal-mass SCs ("equal-mass" case) and, (2) the SC system is composed of SCs with different 
masses ("multi-mass" case). For (1), each SC is assumed to have a mass of 2 x 10^ Mq (or 
My = —10.35 mag) and asc=6.8pc. For (2), the SCs are assumed to have a luminosity (M) 
function consistent with that observed: 

<I>(M) = constant x e-(^-^«)'/2-™', (5) 

where Mo(l^)=— 7.27 mag and o"m=l-25 mag (Harris 1991). The number of SCs (A^sc) in 
a system is a free parameter, ranging from 2 to 20 for the equal-mass case and 2 to 200 
for the multi-mass one. The model with the maximum A^sc of 200, corresponds to the most 
massive VMSCs that have been observed (Drinkwater et al. 2003). In our simulation, the 
masses of the stellar particles in the SCs are assumed to be equal, so that the total number 
in the simulation depends on A'^sc. For example, the particle number in the multi-mass model 
with A^sc = 200 is 214858. We mainly describe the equal-mass 3D model with A^sc = 12 and 
Cv = 0.5 (the "fiducial" model), because this model shows both typical behavior of VMSC 
formation and one of the most interesting results in the present study. Also, we describe 
the (1) scaling relations of the simulated VMSCs and (2) parameter dependences of VMSC 
properties, based on 60 different models. The mass, length, time, and velocity units are 2.0 
X 10^ Mq, 34.0pc, 2.1 X lO^yr, and 15.9kms-^ 



3. Results 

As seen from the fiducial model shown in Figure 1, smaller SCs repeatedly merge with 
one another to form bigger clusters through the process of dynamical collapse of the SC 
system. These bigger clusters then merge to form a single VMSC with an outer diffuse 
stellar envelope, within ~ 10^ yr. The VMSC has an effective radius {Re) of 19.4 pc (2.85 
Re of the progenitor SC) and within 5i?e a mass of 2.1 x 10^ Mq, which corresponds to My 
= —12.9 mag and is 10.2 times more massive than the mass of the original SCs. Figure 2 
shows the structural and kinematical properties of the VMSC. Three different non-spherical 
shapes can be clearly seen in the three projected mass distributions, which suggests that the 
VMSC is a triaxial system. The ellipticity (e), defined as e = 1-b/a, where a and b are the 
long and the short axes, respectively, in the isodensity contour of the projected mass profile, 
is estimated to be 0.14 for the x — y projection, 0.08 for the x — z projection, and 0.27 for the 
y — z projection at the effective radius of the VMSC. Due to efficient conversion of the initial 
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orbital energy of the SCs into internal rotational energy during the SC merging process, the 
final VMSC has a non-negligible amount of rotation that is indicated by moderately high 
Kn/<7o (~ 0.3), where (Tq and V^ are the central velocity dispersion and the maximum line- 
of-sight velocity of the VMSC for each projection. It should be noted here that a fiattened 
triaxial system with non-negligible rotation is remarkably different from typical GCs that 
have no net rotation and quite spherical shapes (mean e = 0.07; e.g.. White & Shawl 1987). 

The simulated VMSCs show interesting correlations between their structural and kine- 
matical parameters (Figure 3). Firstly, more luminous VMSCs have larger central velocity 
dispersions (ctq) and this correlation can be expressed as (Tq oc L°'^^, the slope of which is 
similar to that (0.25) of the Faber-Jackson(1976) relation derived for elliptical galaxies. Sec- 
ondly, more luminous VMSCs have larger effective radii and the correlation can be expressed 
as i?e oc L°'^^, though the dispersion in this relation is moderately large. This can be com- 
pared with the corresponding relation {Re oc L^°^) derived for elliptical galaxies (Kormendy 
1985). Thirdly, L is more strongly correlated with the central surface brightness (ho) than 
the half-light averaged surface brightness (Je). Although a single line is fitted to all the data 
points for the L — Je relation, it is possible that the relation is different over the range —12 
< My < —10 mag and —14 < My < —12 mag (i.e., there is a hint of V-shaped distribution 
in Figure 3). 

Figure 4 shows the comparison between the locations of all of the simulated VMSCs in 
the My- o"o plane and the corresponding observations (Drinkwater et al. 2003). Here only (5) 
UCD points are plotted, since data for NCs is not available. The locations of the simulated 
brighter VMSCs are consistent with the observations, and both the simulated and observed 
data points are closer to the Faber- Jackson relation than to the My — o"o relation of GCs 
(Djorgovski et al. 1997). This implies that the origin of UCDs structural and kinematically 
properties is significantly different to that of GCs, and is closely associated with physics of 
multiple merging of SCs. Thus the present results on the scaling relations of VMSCs clearly 
show that the scaling relations of VMSCs formed from multiple SC merging are significantly 
different both from those of their progenitor SCs (or GCs) and from those of dynamically 
hot early-type galaxies. 

The parameter dependences of structural and kinematical properties of the simulated 
VMSCs can be summarized as follows: Firstly VMSCs are likely to be more flattened (e = 
0.2 ~ 0.3) in the 2D models than in the 3D models for given parameter values of Cy and A'^sc. 
Secondly, triaxial VMSCs in some 2D models show large \4i/o"o (~ 0.4) and figure rotation 
like barred galaxies. Thirdly, the multi-mass 3D models with large Nsc (> 50) show both 
smaller e (i.e., less flattened) and smaller Kii/o'o (i-e., less strongly supported by rotation). 
This is because a larger number of SCs merge with one another from random directions in the 
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multi-mass 3D models. These results suggest that the structural and kinematical properties 
of stellar galactic nuclei (i.e., NCs and UCDs) can differ, depending on the merging histories 
of SCs. 



4. Discussions and conclusions 

We have demonstrated that if VMSCs are formed from the multiple merging of SCs, with 
the observed scaling relations of GCs, the scaling relations of VMSCs are very different from 
those of GCs. Ongoing and future photometric and spectroscopic observations (e.g., HST 
ACS and Keck 10m) of the structural and kinematical properties of VMSCs will therefore be 
able to assess the viability of the "cluster cannibalism" scenario of stellar nucleus formation. 

We have also shown that (1) the intrinsic shapes of VMSCs are more likely to be triaxial, 
and (2) some VMSCs can have rotational kinematics. We thus suggest that further obser- 
vations which provide better statistics on: (l)the e distributions of the projected isophotal 
shapes of VMSCs (which strongly depend on the intrinsic shapes of VMSCs), and (2) the 
locations of VMSCs on the e — Kn/o"o plane (which depend on the internal kinematics of 
VMSCs) will help discriminate between different VMSC formation scenarios. 

Our study also has other important implications. Firstly, the significant amount of 
rotation observed for the metal-poor populations in the most massive Galactic globular 
cluster uj Cen (e.g., Norris et al. 1997), provides evidence that u Cen may have originated 
from the nucleus of a dwarf galaxy. A growing number of observations and theoretical studies 
have recently suggested that oo Cen was previously a nucleus of a ancient nucleated dwarf 
orbiting the young Galaxy (e.g., Hilker & Richtler 2000, 2002; Bekki & Freeman 2003). The 
present study has demonstrated that VMSCs formed from SC merging can have a significant 
amount of rotation (i.e., larger Kn/o"o)) due to the conversion of orbital energy into intrinsic 
rotational energy. Therefore the rotational kinematics of the metal-poor populations of 
UJ Cen reflect the past merging of SCs in the central region of its host galaxy. 

Secondly, we can significantly underestimate the mass-to-light-ratios (M/L) of VMSCs, 
in particular those with large V^/ao, if we estimate the M/L by adopting the commonly used 
formula in which M/L is derived only from the central velocity dispersion and effective radius 
(e.g., Meylan et al. 2001). This is simply because the implicit assumption (in the formula) 
that kinematical energy can be accurately measured by the central velocity dispersion ctq 
alone is not valid for these VMSCs. Accordingly, the real values of M/L for UCDs may be 
even larger than the moderately large values {M/L = 2-4) that are observed (Drinkwater et 
al. 2003), if UCDs indeed have rotational kinematics - something which future spectroscopic 
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observations of UCDs can confirm. 

Thirdly, the observed young stellar populations in NCs (e.g., Boker et al. 2004b) can 
be due to nuclear star formation triggered by dynamical interaction between a NC and its 
surrounding interstellar gas. Recent numerical simulations have suggested that dynamical 
interaction between self-gravitating triaxial systems with figure rotation and the surrounding 
gas can cause rapid gas transfer to the central region and trigger subsequent starbursts there 
(Bekki & Freeman 2002). The present study has shown that some triaxial VMSCs have 
figure rotation (particularly when they are formed from SCs with initial disky distributions). 
Therefore, we suggest that the observed young populations in NCs may be due to past 
dynamical interaction between the triaxial NCs with figure rotation and the surrounding 
gas. 

Although we have not discriminated between UCDs and NCs throughout this paper, 
there exists a possible remarkable difference in the size distribution between the two: The 
observed effective radii range from 15 pc to 30 pc for UCDs and from 2.4 pc to 5pc for NCs. 
This difference can possibly be ascribed to the young stellar populations with compact spatial 
distributions within NCs (i.e., the size distributions for old populations of NCs could not 
be so different from those of UCDs). However, it is equally possible that the difference is 
due to the differences in formation processes between UCDs and NCs. Differences in scaling 
relations between UCDs and NCs, if any, would tell us something about the origin of the 
difference between UCDs and NCs. 

K.B. and W.J.C. acknowledge the financial support of the Australian Research Council 
throughout the course of this work. 
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Fig. 1. — Morphological evolution of twelve equal-mass star clusters (SCs) projected onto 
the X — z plane for the fiducial model (with an initial total mass of 2.4 x 10^ ^o)- The 
length and the time are given in our units (34.0 pc and 2.1 x 10^ yr, respectively). Here the 
time T represents the time that has elapsed since the start of the simulation. 
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Fig. 2. — Final mass distributions {upper) and kinematical properties {lower), projected onto 
the X — y plane {left), the x — z plane {middle), and y — z plane {right) in the fiducial model 
at T = 32.0 in our units (corresponding to 6.7 x 10^ yr). Solid and dotted lines in each lower 
panel represent the radial profile of line-of-sight velocity and that of the velocity dispersion, 
respectively. The length and the velocity are given in our units (34.0 pc and 15.9 km s~^, 
respectively). 
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Fig. 3. — Correlations of structural and kinematical parameters with My (V— band absolute 
magnitude) for the VMSCs in 40 models including equal-mass 2D and 3D models with Cy 
=0 and 0.5. Projected central velocity dispersion (ctq; upper left), half-light-averaged surface 
brightness (Jg; upper right), effective radius {Re] lower left), and central surface brightness 
(Jio; lower right) are plotted against My. Here the central surface brightness Jio is expressed 
as O.IL/tt/Rk)^, where L, Rio are the total luminosity of a VMSC and the radius within 
which 10% of L is included, respectively. The best fit scaling relation for the VMSCs is 
derived for each panel using the least square fitting method and described as a dotted line 
with the derived relation (e.g., ao ex L°-^^). 
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Fig. 4. — Correlations of (Tq with My for the simulated VMSCs {filled squares) and the 
observations {crosses). The results of 60 models, including both equal-mass 2D/3D models 
and multi-mass ones, are shown. Only 5 UCDs with known (Tq (Drinkwater et al. 2003) are 
plotted (no velocity dispersion data for NCs are available). For comparison, the observed 
relations are given by a dashed line for GCs (Djorgovski et al. 1997) and a dotted line for Es 
(Faber & Jackson 1976). 
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